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ABSTRACT 

The effective extinction law for supernovae surrounded by circumstellar dust 
is examined by Monte-Carlo simulations. Grains with light scattering properties 
as for interstellar dust in the Milky- Way (MW) or the Large Magellanic Clouds 
(LMC), but surrounding the explosion site would cause a semi-diffusive propa- 
gation of light up to the edge of the dust shell. Multiple scattering of photons 
predominantly attenuates photons with shorter wavelengths, thus steepening the 
effective extinction law as compared to the case of single scattering in the in- 
terstellar medium. Our simulations yield typical values for the total to selective 
extinction ratio Ry ~ 1-5 — 2.5, as seen in recent studies of Type la supernova 
colors, with further stiffening differential extinction toward the near-UV. 

Subject headings: dust, extinction — supernovae: general 



1. Introduction 

The uncertainties in the brightness corrections of Type la supernovae (SNIa) for color 
excess is among the largest syst ematic uncertaintie s in the use of Type la supernovae to 
measure cosmological distances (INordin et al.l 120081 ). The standard interpretation of color 
excess being due to extinction by interstellar dust in the supernova host galaxy has recently 
been challenged by the empirically deduced color-brightness relation for SNIa. While stud- 
ies of differential extinction of quasars shining through foregrou nd galaxies yield va lues of 
Rv = Av/E{B — V) compatible with the average MW value ( Ostman et al. 2008 1^. the 
SNIa Hubble diag r am scatter is minirn ized for values significantly smaller than Ry = 3.1 
( lAstier et al.ll2006l : iKowalski et alll2008l ). Furthermore, the use of SNIa for cosmology bene- 
fits significantly from the understanding the extinction law for the full optical and near-IR 



shou ld be noted that lo w Ry values to individual QSO systems have been found in e.g. (jOstman et al. 

2008h and (|Wang et al.ll2004h 
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wavelength range. E.g. the wavelength dependen ce of extinction tow ard the Magellanic 
Clouds differ from the Milky- Way extincti on law in (ICardelli et al.lll989[) . even for very sim- 
ilar values of Ry- For quasar sight-lines, lOstman et al.l (120081 ) tested both Milky- Way like 
extinction as well as Small Magellanic Cloud (SMC) extinction law, both giving comparable 
goodn ess of fit. A preference for SMC dust for extinction of AGNs has been suggested bylLj 

feoorh . 



Recently, the detection of circumstellar (CS) matter in the local environment surround- 
ing the Type la supernova SN2006X in the nearby galaxy MlOO has been reported by 



Patat et al.l (120071 ). A shell within a few 10 cm (~ 0.01 pc) of the center of the explo- 



sion has been suggested to explain the time-variable Na I D lines in the SN spectrum. 



Wang et al.l (l2008al ) report Ry = 1.48 ± 0.06 and E( B - V) = 1 . 42 ± 0. 04 mag for 



SN2006X and a light echo in the hghtcurve was found by IWang et al.l (l2008bl ) consistent 
with dust illuminated at a distance of 27-170 pc from the site of the explosion. Even if 
the local environment around this supernova may not be very common among SNIa, similar 
values for the total to selective e xtinction ratio have be en reported for several SNIa with 
good w avelength coverage. E. g . iKrisciunas et al.l (120071 ) found Ry = 1-55 ± 0.08 for SN 



1999cl; Elias-Rosa et al. 



J2006I . l2008h report Ry = 1.80 ± 0.19 and Ry = 1.59 ± 0.07 for 
SN2003cg and SN 2002cv respectively. Furthermore, a s tatist ical study of optical colors 
of a sample including 80 near-by SNIa, iNobili fc Goobarl (120081 ) found an average value of 
Ry = 1.75 ± 0.27 for SNIa with E{B — V) <0.7, and even lower for a subsample of low- 
reddening SNIa. 

Next, we examine the possibility that low values of Ry stem from the semi-diffusive 
propagation of photons in the neighborhood of the site of the supernova explosion. 



2. Quasi-diffusive light propagation around the supernova 

Photon propagation around a mediu m of scatterers can be described by a quasi random- 



walk picture. The reader is referred to (IChandrasekharl Il943l ) for a beautiful introduction 



to this subject. Lets consider a localized distribution of dust particles within a distance 
Res from the explosion site, negligibly small compared to the distance to the observer, d, 
i.e. RsN < Res ^ d, where Rsn corresponds to the radius from where the SN radiation 
emerges. The trajectory of a photon will be straight until it hits a dust particle at which point 
the photon can either be scattered or absorbed. If the photon is scattered in a new direction, 
it follows a straight path until the next encounter, and so on until r > Res- The mean free 
path between interactions, Ag//, is thus determined by the number density of scatterers, n, 
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and their effective cross-section for scattering and absorption of liglit, cTg// = as + (Ja- 

1 



A 



eff 



n ■ aeff 



For Ae// ^ Res J the single scattering approximation is valid and the light beam reaching 
the observer is attenuated as e~^, with r = Rcs/Kff- This is the case applicable for 
extinction by dust in the interstellar medium. 

For Ae// ^ Res, corresponding to a (local) high number density of scatterers, the situa- 
tion is different. If the absorption probability is much lower than for scattering, photon prop- 
agation is diffusiv e and the average pro perties can described analytically using the formulas 
for random-walk (lAskebjer et al.lll997l ). The case we are considering here is for Ae// ~ Res, 
where the scattering cross-section exceeds the absorption cross-section , > cTq, in a wave- 
length dependent manner. In particular, we examine the cases where the light scattering 
properties of dust particles in the CS matter are si milar to what has been modeled for 
interstellar dust grains the Milky- Way (IDraind l2003l ) or the LMC ( IWeingartner fc Draine 
200lh . 



Table d] shows the the wavelength dependent albedo factor (= as/{as + cXa)) and the 
average of the cosine of the scattering angle for interactions between light and dust particles. 
Also tabulated is the absorption cross-section divided by dust mass. Note that the Milky- 
Way parameters corr espond t o a d ust size distribution matching Ry = 3.1 for dimming of 
stars in the Galaxy. IWang (120051 ) considered the impact of circumstellar dust upon the 
measured value of Ry, but only the extreme case where all scattered photons may reach the 
observer. That assumption overlooks an important aspect of the problem: while the bluer 
photons scatter more, they also are more likely to be absorbed. This leads to a steeper 
wavelength dependence of the effective extinction law, possibly explaining the unusual total 
to selective extinction ratios found in studies of SNIa. 



3. Monte-Carlo simulation of light propagation around supernova 

In order to estimate the net effect of scattering and absorption on the light reaching 
the outer edge of a shell of circumstellar dust around the SN site. Res, a Monte-Carlo 
simulation was performed. Photons with energies corresponding to the central wavelengths 
of the UBVRIJHK photometric system were generated and subsequently followed as they 
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propagate in the dusty mediumj^ A uniform distribution of scatterers within a sphere of 
radius Res is used in the calculations. Our treatment is rather insensitive of the physical size 
of Res since what governs the differences in path-lengths of photons at different wavelengths 
are the optical depths, = Rcsl\ and = Rcs/^a- Thus, for a fixed color excess, a larger 
Res can be compensated by a lower number density, n, thus keeping Ts and Tq unchanged. 
Figure [1] shows the wavelength dependence of and in our calculation for MW and LMC 
dust types for a reddening at r = Rqs of E{B — V) = 0.1. The key observation is that 
at short wavelengths, both and increase, while dominates. Since they scatter more, 
photons with short wavelength leaving the CS region to eventually reach the observer must 
propagate a larger path-length than photons at longer wavelengths. However, since also 
increases at shorter wavelengths, the relative attenuation of bluer photons is enhanced, thus 
generating an effective extinction law with steeper wavelength dependence, i.e. lower Ry- 



SMC-type dust, as parameterized in (IWeingartner fc Drainell2001l) . has a smaller albedo 
in the optical to near-IR region (see Fig. 23 in IWeingartner &: Draind (120011 )). leading to a 
larger value of compared to LMC or MW dust. Thus, SMC dust is not suitable to explain 
the anomalous SN colors with the scenario presented in this work. 

For each photon starting at Rsn ^ Res, straight path- lengths between interactions 
{Ls,La) were generated from an exponential distributions A^^ exp {—Lg/ Xs) and A^^ exp {—La/ Xa), 
where the mean free path for scattering and absorption, Xg = {n ■ o"s)~^ and Xa = {n ■ 0"^)^^ 
were calculated from the wavelength dependent parameters in Table [H The scattering an- 
gle for each int eraction was generated in the simulation following the Henyey-Greenstein 
approximation (IHenyey fc Greensteirull94ll ): 



da 



d(cos(^)) {I + g2 - 2g cos{e)) 



3 ' 
2 



(2) 



where g =< cos{6) > (also listed in Table [1]). The density of dust in the shell was varied 
to cover a wide range of reddening at r = Res, < E{B — V) < 0.4. The probability for 
photons to reach Res without absorption is calculated by repeating the ray-tracing Monte- 
Carlo 5 ■ 10^ times for each wavelength. 



^In this study we neglect the effect noted by IWang (|2005| ) where the effective wavelength in each filter is 
changing with time following the color evolution of SNIa 
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Low Rv from simulations 



Figure [2] shows that the attenuation of light after the circumstellar shell yields low values 
of Rv = Ay/ E{B — V), as anticipated. In particular, the simulation s using LMC dust , resul t 
in an effective Ry = 1.65, compatible with the b est fit result of (INobili &: Goobarl |2008| ) . 
We note, however, that the CCM extinction law (ICardelli et al.lll989l ) does not accurately 
reproduce the wavelength dependence of photons attenuation over the entire optical range, 
as shown in Figure [^. The deviations are sig nificant for the U-B color, in agreement with 
the observations reported in (iGuy et al.ll2005l ). 



Instead, a power- law relation gives a good fit to the simulations (see Figure Hj): 

—— = 1 - a + a — 
Ay \Xy 



(3) 



where Xy = 0.55 /xm is the central wavelength of the V-band filter. Since the attenu- 
ation at B-band (central wavelength A^) is given by {Ry + 1)E{B — V), the parameters in 
Equation ([3]) are related to Ry as 



R 



V 



1 



-1 



a(0.8P- 1) 



(4) 



It is encouraging for the use of Type la SNe for precision cosmology that the reddening 
corrections may take a simple and general analytical form, if the model is confirmed. Further 
studies of SN colors, also including near-IR data, should be used to test this relation. 



5. Summary and conclusions 



Simple simulation s show that circu r nstellar material, d etected in at least one Type la 
supernova, SN2006X (jPatat et al.l 120071 : IWang et al.l l2008bl ). could potentially explain the 
empirically determined extinction law for low redshift SNIa, especially if the circurnstella r 
material resembles LMC dust grains. Adopting the CS shell size of iPatat et al.l (120071 ). 
10^^ cm, we find that for r ~ 1, the required mass in dust around the supernova 



Res ~ 
is Mdust 
cross-section from Table [1] 



^T^RcsKcTa/mdust) 



10 ^Mq when inserting typical values of the absorption 



•^The validity of the simulation code was tested by accurately reproducing the standard extinction laws 
when the diffusive process was switched off. 
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A simple power-law expression is found to fit very well the effective extinction law for 
dust in the CS environment of the supernova produced by Monte-Carlo simulations. 

Depending on the thickness of the CS shell, shifts in the time of lightcurve maximum 
may be expected for different bands since the amount of quasi random-walk will differ. In 
particular, photons in redder bands will suffer less scattering and thereby less time delay. 
This effect should correlate with the measured reddening, E{B — V). The assumption of 
a uniform density is not expected to be critical for the results at first order. However, a 
second order effect may be expected since a large scale of Res would result in a longer time 
for photons being "trapped" in the scattering sphere. As the the intrinsic colors of Type la 



change on a time scale of days (INobili fc Goobarll2008l ). time delays of photons of that time 



scale would affect the measured colors as a function of time. 

In a forthcoming paper, potential direct observables from interaction between photons 
and dust will be investigated. Also, the sensitivity of the effective extinction law to the dust 
grain sizes and density profile in CS medium and the combination of both scattering in the 
circumstellar material and the interstellar medium needs to be further investigated. If the 
presence of CS material is indeed the source of the color-brightness relation found in SNIa, 
the case for restframe near-IR observations is further strengthened: the peak magnitude 
corrections, and their model dependence, are smaller than at optical wavelengths. 



I am grateful to Vallery Stanishev for discussions triggering this work and to the anony- 
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manuscript. 
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Fig. 1. — Optical depth for scattering (r^ = Rcs/K) and absorption (tq = Rcs/^a) as a 
function of wavelength from the Monte-Carlo simulation of photon propagation in a dusty 
circumstellar material around t he SN for E( B — V) = 0.1 at r = Rns- Scattering proper 



ties for average MW dust from (jDraindl2003l ) and LMC dust (jWeingartner &: Draindl200ll ) 
tabulated in Table [U were used. The bluer photons scatter more and since they also have 
a larger probability for absorption, fewer make it to the edge of the dusty shell, Res- As a 
result, a steeper wavelength dependence for attenuation is to be expected. 
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Fig. 2. — Rv = Av/E{B — V) vs E{B — V) resulting from the Monte-Carlo simulation 
of photon propagation in a dusty circumstellar material around the explosion site. The 
"wiggles" in the curve are compatible with the statistical un certainty of t he Monte-Carlo 
simulation. Scattering prope rties for average MW dust from (jDrainell2003l ) and LMC dust 
( IWeingartner &: Drainell200ll ) tabulated in Table [H were used. The average total to selective 
extinction ratios found were Ry = 2.56 for Milky- Way type dust and Ry = 1-65 for dust 
compatible with properties in the Large Magellanic Clouds. 
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Fig. 3. — The extinction law, Ax/ Ay, derived from the Monte-Carlo simulat ion of circum- 



stella r dust are compared with the parameterized Milky- Wa y extinction la w (jCardelli et al. 



19891). scattering properties for average MW dust from (IDraind 120031 ) and LMC dust 



Weingartner fc Draind l200ll ) tabulated in Table [T] were used. The two cases considered 
yield a poor match to the standard Ry = 3.1 case. Even when the Ry is adjusted to fit 
the Av/E{B — V), differences are found in the wavelength dependence, especially at shorter 
wavelengths. 
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Fig. 4. — Empiri cal functiona l fit to the simula ted data, scattering properti es for average 
MW dust from (jDraind |2003| ) and LMC dust (IWeingartner h Draind I2OOII ) tabulated in 
Table [T] were used. 
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Table 1. Scattering parameters for optical and near-IR photons from 



(jWeingartner fc Drainell200ll : lDraindl2003l ) corresponding to interstellar extinction in the 
Large Magellanic Clouds (LMC) and Milky- Way (MW) with i?y = 3.1. 



wavelength 


albedo 


g=< cos(6) > 






filter 


(fira) 


MW 


LMC 


MW 


LMC 


MW 


LMC 




0.36 


0.6203 


0.6532 


0.5695 


0.6072 


1.595E+04 


1.162E+04 


U 


0.44 


0.6529 


0.7159 


0.5654 


0.6153 


1.191E+04 


7.542E+03 


B 


0.55 


0.6735 


0.7631 


0.5382 


0.6059 


8.551E+03 


4.666E+03 


V 


0.65 


0.6745 


0.7785 


0.4995 


0.5815 


6.722E+03 


3.319E+03 


R 


0.80 


0.6576 


0.7686 


0.4381 


0.5289 


5.170E+03 


2.416E+03 


I 


1.22 


0.5846 


0.7029 


0.2893 


0.3892 


3.224E+03 


1.550E+03 


J 


1.63 


0.5200 


0.6472 


0.2086 


0.2911 


2.351E+03 


1.237E+03 


H 


2.19 


0.4391 


0.5599 


0.1310 


0.1792 


1.670E+03 


1.041E+03 


K 



